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Abstract 

It is shown that there exist families of asymptotically flat solutions of the Einstein 
equations coupled to the Vlasov equation describing a collisionless gas which have a New- 
tonian limit. These are sufficiently general to confirm that for this matter model as many 
families of this type exist as would be expected on the basis of physical intuition. A central 
role in the proof is played by energy estimates in unweighted Sobolev spaces for a wave 
equation satisfied by the second fundamental form of a maximal foliation. 
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1. Introduction 

It is a well known empirical fact that in many situations a general relativistic descrip- 
tion of the motion of self-gravitating matter can be replaced to a good approximation by 
a non-relativistic, Newtonian one. In the usual formulation of Newtonian gravity the in- 
teraction is described by a single scalar function, the Newtonian potential. The relation of 
this to general relativity, where the fundamental object is a Lorentz metric, is obscure. The 
basic idea required to understand this relation mathematically was provided by Cart an [4]. 
He showed that Newtonian theory can be formulated in such a way that the basic object 
is an affine connection whose non-zero components are components of the gradient of the 
Newtonian potential. The role of the potential itself is then merely that of providing a 
convenient representation of this connection in certain coordinate systems. It was realised 
by Friedrichs[12] that the natural way to connect the two theories is to require that the 
Levi-Civita connection of the spacetime metric go over in the limit as the speed of light 
c goes to infinity into the connection defined by Cartan. Since then many authors have 
extended this work on the relations between the equations of the two theories and the 
physical interpretations of their solutions. This knowledge has been systematised in the 
frame theory of Ehlers (see [11] and references therein). What has been achieved is to 
set up a precise definition of the Newtonian limit of general relativity which encodes that 
which is desirable on physical grounds. The major open question is whether this definition 
is compatible with the Einstein equations in the sense that there exists a sufficiently large 
class of solutions which satisfy all the axioms. The purpose of this paper is to answer this 
question in the affirmative. 

The case of the Newtonian limit which is of most physical interest is that of an isolated 
system. This is expressed mathematically by restricting attention to asymptotically flat 
solutions of the Einstein equations. (The case of cosmological solutions, which is also of 
considerable interest, will not be treated here.) It is of prime importance to have results 
which do not only handle the vacuum Einstein equations since in that case only the trivial 
Newtonian solution (i.e. empty space) could be expected to limit of singularity 

free asymptotically flat spacetimes. It is at this point that the first serious difficulty is 
encountered. The most obvious type of matter to take would be one or more bodies made 
of fluid or an elastic solid. However only very limited results exist on the initial value 
problem for self-gravitating bodies of this type [20]. It is for this reason that the matter 
model chosen here is a collisionless gas described by the Vlasov equation. It is known that 
the local in time initial value problem for the Vlasov-Einstein system is well posed for a 
class of initial data which allows spatially localised matter distributions [5]. A large part 
of what follows does not crucially depend on any property of the particular matter model 
chosen beyond the fact that the local in time initial value problem is well posed. There 
is, however, one step where a special property of the Vlasov equation is used, namely 
in the last paragraph of section 3. In order to generalise the results of this paper to a 
different matter model it would be essential to find a replacement for the argument of that 
paragraph or to modify the structure of the main proof significantly to avoid the need for 
that argument. 

The next difficulty which hampers the development of rigorous theorems on the New- 
tonian limit is that this limit is singular in the sense that the Einstein equations, which are 
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essentially hyperbolic, go over into the Poisson equation, which is elliptic. The hyperbolic 
nature of the Einstein equations is related to the propagation of gravitational waves. There 
is one special case, namely the case of spherical symmetry, where gravitational radiation is 
absent. This leads to a simplification of the problem and the Newtonian limit of spherically 
symmetric asymptotically flat solutions of the Vlasov-Einstein system was handled in [19]. 

When spherical symmetry is not assumed the problem of the singular limit has to 
be faced and to see which way to go it is useful to consider a simpler analogue of the 
Vlasov-Einstein system where a limiting situation occurs which is rather similar. This is 
the Vlasov-Maxwell system whose quasi-static limit has been considered in [1] , [8] and [22] . 
Of these papers the one which is of most relevance here is that of Degond[8]. He treats the 
limit using the fact that in the energy estimates for the Maxwell equations the terms in the 
equations which blow up as c — > oo make no contribution. The solutions discussed belong 
to a Sobolev space on each slice of constant time. For the Einstein equations this does not 
hold. An asymptotically flat metric falls off only as r _1 as r — > oo on a spacelike slice and 
the positive mass theorem implies that any attempt to impose faster fall-off excludes all but 
the trivial solution. Thus the metric does not belong to a Sobolev space. The usual way to 
get around this is to replace the ordinary Sobolev space by a weighted one. Unfortunately 
it is easily seen that such a replacement destroys the property used by Degond that singular 
terms drop out of the energy estimates. In the following this difficulty is circumvented by 
using a formulation of the Einstein equations due to Christodoulou and Klainerman[6]. 
There the only object which is determined by solving a non-trivial hyperbolic equation is 
the second fundamental form, which does lie in an ordinary Sobolev space. The r _1 part 
of the metric is generated by solving an elliptic equation. 

Now that the strategy has been outlined, the main theorem will be stated. The 
notation is as follows: g ab is the induced metric on the leaves of a maximal foliation of 
spacetime, k ab is the second fundamental form of this foliation, is the lapse function, 

are the Christoffel symbols of the spacetime metric and / is the phase space density 
of particles. The spacetime metric is of the form 

ds 2 = -cj) 2 dt 2 + g ab dx a dx b . (1.1) 

The notions of regular initial data and regular solutions appearing in the statement will 
be defined in Sect. 5. The exact interpretation of the order symbols used will be given at 
the end of Sect. 6. The parameter A corresponds physically to c~ 2 . 

Theorem 1.1 Let (<7° ft (A) , k® b (\) , / '(A)) be a parameter- dependent initial data set for 

the Vlasov-Einstein system which is regular of order s for some s > 6 and satisfies the 

constraints and the maximal slicing condition. Suppose that as A — > 0; 

(i) g° ab (X) = X8 ab + 0(A 3 / 2 ) 

(it) k° ab (X) = 0(A 3 / 2 ) 

(in) d t k° ab (X) = 0(A 3 / 2 ) 

(iv)f(X) = f N + o(l) 

for some fy. Then a solution (g ab (X),k ab (X),(f)(X), f(X)) of the Vlasov-Einstein system, 
which is regular of order s and induces the given initial data on the hypersurface t = 0, 
exists on a X-independent time interval [0, T) and has the properties: 
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(i) 9ab(X) = A5 Q 6 + o(A) 
(%%) k ab (X) = o(A) 

(ni) 0(A) = 1 — £7A + o(A) /or some £/ 

^;rg (A) = -5 ab V6t/ + o(i) 

(v,) a// ot/ier components are o(l) 
(vi)f(X) = f N + o(l) /or some /at 

Moreover /jv and U solve the Vlasov-Poisson system with initial datum f^. 

In assumption (iii) in the hypotheses of this theorem the time derivative is to be calculated 
using the Einstein evolution equation (2.4) and the function <f> in that equation is to be 
got by solving the lapse equation (2.12). This assumption may seem unnatural but it is 
essential. Its significance will be discussed further in Sect. 7. 

It is appropriate at this point to mention some recent work related to the present paper. 
Fritelli and Reula[13] have suggested an interesting approach to proving convergence of 
solutions of the Einstein equations in the Newtonian limit on a spatially bounded region. 
Lottermoser[15] has proved the existence of rather general families of solutions of the 
Einstein constraint equations having a Newtonian limit. It was necessary to prove some 
new results on existence of families of solutions of the constraints in the present paper since 
Lottermoser's method is not suitable for producing solutions which satisfy a prescribed 
gauge condition (e.g. the maximal slicing condition used in the following). It is also of 
interest that the solutions whose existence is demonstrated in the present paper include 
ones which do not belong to the class produced in [15]. This is because the basic object Z ab 
used there and supposed to behave regularly as A — * diverges in general in the Newtonian 
limit for the data constructed here. 

The paper is organised as follows. In Sect. 2 the form of the Einstein equations used 
in [6] is discussed and it is shown how the parameter A can conveniently be introduced 
into it. In the third and fourth sections estimates are derived for the Vlasov and Einstein 
equations respectively. In Sect. 5 these are used to prove the local existence of a solution 
on a A-independent interval and the remainder of Theorem 1.1 is proved in Sect. 6. The 
existence of a large class of regular initial data is demonstrated in Sect. 7. Two appendices 
are concerned with some elliptic theory and estimates for modified Sobolev spaces which 
are needed in the body of the paper. 



2. Derivation of the reduced equations 

Consider first the 3+1 form of the Einstein equations with zero shift. The constraints 

are 

R- \k\ 2 + (trk) 2 = 16n<p- 2 T 00 , (2.1) 

V a k ab - V b tvk = Sn^Tot, (2.2) 

and the evolution equations are 

d t g a b = -2#ab, (2.3) 
d t k a b = -V a V 6 + 4>{R a b + trkk ab - 2k ac kl 

- 8ivT ab - 47r0- 2 Too#a6 + 4irTg ab ) . (2.4) 
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Here R ab is the Ricci tensor of g ab . The objects T 00: T 0a and T ab are components of the 
matter tensor and T is obtained by taking the trace of T ab with the metric g ab . Define 

A = trfc, (2.5) 
B = R-\k\ 2 , (2.6) 

C a = V b k ab -^V a (trk), (2.7) 

D ab = (p^dtkab + <T 1 V a V 6 </> - R ab + 2k ac k c b . (2.8) 

It will also be useful to have the modified quantities 

B = B -16n<p- 2 T 00 , (2.9) 
C a = C a + 8Tr<f)- 1 T 0a , (2.10) 
Dab = D ab + 8nT ab + 4n(<p- 2 T 00 - T)g ab . (2.11) 

In the following the maximal slicing condition A = will be used. Under that assumption 
equations (2.1), (2.2) and (2.4) are equivalent to B = 0, C a = and D ab = respectively. 
If A = then the lapse equation 

A0= (\k\ 2 + 4ir(j)- 2 Too + 47rT)(j) (2.12) 

is satisfied. Following [6] it can be shown that (2.1)-(2.4) together with A = imply a 
wave equation for k ab . 

-(^dtfkab + Ak ab = N ab + r ab , (2.13) 

where 

~1 



T ab = 8lT 



<t>- 3 T 00 gab - ^(VaTob + V b T Qa ) + 0" 1 T a6 



-]^- l fg ab + cp- 2 (To a V b( p + T 06 V a </>) - ((p- 2 T 00 - T)k ab - <\>-^T m g ab 



(2.14) 



N ab = L ab - H ab , (2.15) 
<p 2 L ab = V a V 6 </> - (T^VaV^ - f c ab V c( p + 2cj>(k c a d t k bc + k c b d t k ac ) 

+ 4<p 2 k ac k cd k bd , (2.16) 
(j)H ab = <j>I ab + V c 0(2V c /c a6 - V a k bc - V b k ac ) - V a 4>V c k bc - V b 4>V c k ac 

- V c V b (pk c a - V C V #6 + A(j>kab (2.17) 

<pi ab = -3(R ac k c b + R bc k c a ) + 2g ab (k cd R cd ) + k ab R. (2.18) 

The reduced system of Einstein equations which will be used consists essentially of (2.3), 
(2.12) and (2.13). Unfortunately the occurrence of the Ricci tensor in I ab causes trouble 
with the existence theory and so it will be replaced by using the relation 

R ab = ^dtkab + 0- 1 V a V 6 + 2k ac k c b + 8nT ab + 4n(<p- 2 T 00 - T)g ab (2.19) 
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which is equivalent to D ab = 0. The following lemma was proved in the vacuum case in 
[6]. 

Lemma 2.1. Let (g a b, ^at> 7 4>) be a solution of (2.3), (2.12) and (2.13) (with (2.19) having 
been substituted into (2.13) to eliminate R a b)- Then if the data (g a b(^)jk a b(0),d t k a b(0)) 
are such that A, B, C a and D ab vanish for t = and if V Q T Q ^ = then A, B, C a and 
D a b vanish everywhere so that (g a b, k a b, 4>) defines a solution of the Einstein equations for 
which the hyper surf aces t= const, are maximal. 

Proof This will only be sketched since it is very similar to the vacuum case. Equations 
(2.3), (2.12) and (2.13) imply (with F = B + trD): 



<jr x d t A = F:=F- 4tkI)- 2 T o - 4ttT, (2.20) 
4T x d t F = AA - 40 -1 V a 0Ca + (P'^R - A^)A 

- 87r[-20- 1 V a T Oa + 2<p- 2 V a <j)T 0a - Ijl^f - 2k ab T ab 

- 30- 4 0T oo + (-<T 2 7oo + T)A + 3/2^- 3 T o] (2.21) 
^dtCa = V h D ab - l/2V a (trD) + ^V^Dab - 1/20" 1 V a <f>F 

- ^AV^Kb - V b Ak ab + Sn^ToaA (2.22) 
^dtDab = V a C b + V b C a (2.23) 

In the derivation of (2.22) the vanishing of the 4-dimensional divergence of T a/3 has been 
used. The latter can also be used to simplify the equation obtained from (2.21) by substi- 
tuting for F and C a in terms of F and C a . Equations (2.20)-(2.23) imply wave equations 
of the form 

- (4>~ 1 d t ) 2 A + AA = M, (2.24) 

— (4>~ 1 d t ) 2 F + AF = M', (2.25) 
-(<p- 1 d t ) 2 C a + AC a = M^ (2.26) 

where M, M' and M" are linear expressions in the quantities A, V 'A, V 2 A, F, VF, 
4>~ l dtF, C, VC, 4>~ 1 d t C : -D, V-D. The proof of the lemma can now be completed by 
deriving an inequality of the form 

S(t) <K I £{s)ds (2.27) 
Jo 

where 

£= f (\A\ 2 + \VA\ 2 + \F\ 2 + \VF\ 2 + U- 1 d t F\ 2 

Jr3 (2.28) 

+ |Cf + |VCf + l^dtCl 2 + \D\ 2 )dV g , 

and dV g is the volume element on R 3 associated with the metric g ab . This would be 
a straightforward consequence of the usual energy inequalities for wave equations if it 
were not for the occurrence of the quantities VD and V 2 A on the right hand side of the 
equations. Note that these do not appear in the definition of £. In fact VD only occurs 
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in M" and V 2 A only in M' and the problem can be overcome as follows. One of the 
terms which needs to be estimated is schematically of the form J Q [j- R3 (dtC'VD)(s)]ds. 
Integrating by parts in time converts this into the sum of a spacetime integral and a 
boundary contribution on the hypersurface labelled by t. The spacetime integral can now 
be handled by a partial integration in space. A partial integration in space should also be 
applied to the boundary term. It is then schematically of the form J R3 (VCD). This can 
be estimated by an expression of the form 

where if is a constant and rj may be chosen to be any positive real number. Choosing it 
so that Ki] < 1/2 we can absorb the first term into S. To handle the second term express 
it as the integral from to t of its time derivative and use equation (2.23). The term 
containing V 2 A can be estimated in an analogous way. Applying Gronwall's inequality to 
(2.27) now completes the proof. □ 

Lemma 2.1 shows that providing we are dealing with a matter model which guarantees 
that V a T a/3 = (and this is in particular true of matter described by the Vlasov equation) 
then solving the reduced system consisting of (2.3), (2.12) and (2.13) suffices to solve the 
Einstein equations. In this paper the unknowns (g a b, k a b, <fi) are time-dependent objects on 
R 3 . Normally, if asymptotically flat situations are to be studied, the boundary conditions 
9 ah — > 8 a b and — > 1 as \x\ — > oo will be imposed. In order to study the Newtonian limit 
the first of these will be replaced by g a b — > X5 a b where the parameter A corresponds to 
c~ 2 . The Newtonian limit then corresponds to the limit A — > and if this is to be regular 
it will be the case that g a b = 0(A) and k a b = 0(A) as A — > 0. In the terminology of the 
frame theory (see [11]) g a b is part of the temporal metric. In view of this dependence on 
A it is convenient to use the variables "fab = X~ 1 g a b and K a b = X~ x k a b- Then -f a b satisfies 
the standard boundary condition that ^ a b — ■> 5 a b as \x\ — > oo. The basic equations are 

Qtlab = -2<t>Kab, (2.29) 

A 7 = X[\k\ 2 + 4n(j)- 2 Too + 4nT](j) (2.30) 
- (^dtfnab + A- 1 A 7 « a6 = \-\N ab + r ab ). (2.31) 

In equation (2.30) the norm of K a b is defined by 7 a & and not by g a b- 

3. Estimates for the Vlasov equation 

For a discussion of the definition of the Vlasov equation in general relativity see [18]. 
Recall that the phase space density, which is the unknown in this equation, is a real- valued 
function on the mass shell P i.e. the submanifold of the tangent bundle of spacetime defined 
by the conditions g a f3P a p^ = —1 and p° > 0. The manifold P can be coordinatized by the 
spacetime coordinates x a together with the spatial components p a of the momentum. In 
the situation considered in this paper it is identified in this way with R 6 x [0, T). Using the 
3+1 formalism introduced in the previous section the explicit form of the Vlasov equation 
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is 



df 
dt 




(3.1) 



n 5/ 



- <9p a 



where |p| 2 = "f a bP a p b - Attention will be confined to initial data for / which are compactly 
supported in R 6 . Also A will be restricted to belong to the interval (0, A ] for some A > 0. 
No loss of generality results since it is only the limiting behaviour as A — > which is of 
interest here. 

In this section estimates for the solution of (3.1) will be obtained in a fixed background 
geometry. It will be supposed that the quantities T^ c , K a b, 4> and A _1 V a ^> are continuous 
and bounded together with their first derivatives with respect to the spatial coordinates 
x a . Let C\ denote a common bound for these. It will furthermore be assumed that there 
exists a positive constant A such that <p _1 < A and A~ 1 5 a b < 7 a & < A5 a b. It follows that a 
similar estimate holds for "f ab . Equation (3.1) says that / is constant along characteristics 
(in the present context these are the lifts of geodesies to the mass shell) and the assumptions 
on the geometry are enough to guarantee the existence of these characteristics. Let R(t) 
and P(t) denote the maximum values of \x\ and \p\ respectively contained in the support 
of / at time t. Then (3.1) implies estimates of the form P(t) < P(0)(1 + Ct)e and 
R(t) < R(0) + f P(s)ds where the constant C depends only on C\ and A. Next the 
Sobolev norms of / will be estimated under the additional assumptions that 7 a & belongs to 
L°°([0, T), K S (R 3 )) and that K ab and A _1 V a belong to L°°([0, T), H S (R 3 )). Here and in 
the following H s (R, n ) denotes the standard Sobolev space of order s and || \\h<> is the norm 
on that space. The L p norm is denoted by || || p . The space K s (R 3 ) is defined to consist of 
those functions / inL°°(R 3 ) with V/ G H S ~ 1 (R 3 ) with norm \\f\\ K . = ||/||oo + ||V/|| H -i. 
Consider now the norm of f(t) in L 2 (R 6 ). Liouville's theorem implies that the L 2 norm 
of /(£) with respect to the geometrically natural volume form is constant. On the other 
hand, under the assumptions already made on the geometry, this volume form defines an 
equivalent L 2 norm to that of the standard volume form on R 6 . Hence ||/(£)||2 < C for a 
constant C only depending on A and C\. To estimate ||/(t)||i7 a for s > a method used in 
[8] will be adopted. If the Vlasov equation (3.1) is written schematically as Xf = then 
the derivative D s f of order s satisfies an equation of the form X(D s f) = Q s for a certain 
source term Q s . At this stage it is necessary to use the assumptions that 7 a & belongs to 
L°°([0, T), K S (K 3 )) and that K ab and A _1 V a belong to L°°([0, T), H S (H 3 )). Let C 2 be a 
bound for their norms in these spaces. Note that if s > 2 the assumptions made up to now 
imply that -f ab E L°°([0, T), K S (R 3 )) since K s (Tt 3 ) is then a Banach algebra. To estimate 
Q s the facts will be used that (see [14], [16] or [23]): 



for k < s. These hold provided the norms on the right hand side of the inequalities exist. 
They are valid for functions defined on the whole of R n or on a bounded domain. In this 



D k 



\\D k {gh)\\ 2 <C{\\ 

9 II oo 

(gh) - gD k h\\ 2 < C(\\Dg\U\h\\ Hs -i + ||</|| ff .||fc||oo) 



(3.2) 
(3.3) 



8 



section the latter case is the one which is relevant, the domain in question being an open 
subset of R 6 which contains the intersection of the support of / with each hypersurface of 
constant time. There results the estimate 

\\Qs-ih<C\\f\\ H .-i fors>6. (3.4) 

Here C depends only on A, C\ and C2 and the Sobolev embedding theorem in R 6 has 
been used. Combining (3.4), the equation X(D S f) = Q s and the inequality 

II/IIh.<c(||/|| 2 + ||dvI|2) (3.5) 

shows that an integral inequality of the following form holds for s > 6. 

\\f(t)\\ 2 H s-i < 11/(0)11^.-1 +C f \\f{t')\\ 2 Ha -,dt\ (3.6) 

Jo 

where C only depends on A, C\ and C^- In fact C\ can be estimated in terms of C<i and 
||</>| |oo using the Sobolev embedding theorem. 
The matter tensor is defined by 

T^ = - J /2/yV| 7 | 1/2 /Po dp 1 dp 2 dp\ (3.7) 

where 7 denotes the determinant of 7^. This is essentially the energy- momentum tensor. 
It is referred to here as the matter tensor since it has been normalised so that in the 
Newtonian limit T 00 becomes the mass density rather than the energy density. The H s ~ x 
norm of the integrand in (3.7) can be estimated by a constant (depending only on A, C\ 
and C2) times the H s ~ l norm of /. Also for all functions F(x,p) whose supports are 
contained in a given compact set an inequality of the following form holds. 

\\F\\ Hs -i <C||F|| H8 -i, (3.8) 

where F(x) = J F(x,p)dp. Combining this information with (3.6) gives 

w^mi-* < cdimwi^ + f \\f{t')\\ 2 HS -,dt'). (3.9) 

jo 

Suppose next that the time derivatives of 7^, K a b and <fi satisfy estimates of the same kind 
as already assumed for the quantities themselves. Then an argument analogous to that 
above leads to estimates of the form 

\\dtf(t)\\h^ < 11^/(0)11^ + C f \\f{t')\\ 2 Ha -, + \\d t f{t')\\l a . 2 dt\ (3.10) 

Jo 

\\f^(t)\\ 2 H s- 2 < C(\\d t f(0)\\ 2 Hs-* + f 11/(011^- + \\dtf(t')\\ 2 Hs -*dt'), (3.11) 

Jo 
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where the constant C now depends on A, Ci, C2, ||0||<x>> the norm of the time derivative 
of in the space L°°([0, T), i7 s_1 (R 3 )) and the norms of the time derivatives of K ab and 
A _1 V0 in L°°([0, T), if s (R 3 )). 

The derivative D s ~ 1 (d t f) still needs to be estimated. To do this, note that the Vlasov 
equation says that the solution is constant along characteristics. Thus the value of / at 
any point is equal to the value of f° at the point where the characteristic through that 
point intersects the initial hypersurface. Let (X(s, t, x, p), V(s, t, x, p)) be the characteristic 
satisfying X(t,t,x,p) = (x,p) and V(t, t, x,p) = (x,p). Then 

/(*, x, p) = f°(X(0, t, x,p),V(0, t, x,p)). (3.12) 

Hence 



d t f(t, x, p) = D x f°(X(0, t, x, p) , V(0, t, x, p))dX/ds(t, x, p) 

+ D p f°(X{0, t, x, p), V(0, t, x, p))dV/ds(t, x, p) 



(3.13) 



Now the quantities dX/ds and dV/ds are bounded in if s_1 (R 6 ). Assuming that f° is in 
H S (R 6 ) we see that in order to estimate the H s ~ x norm of dtf it suffices to know the 
following two things. Firstly, the mapping taking (t,x,p) to X(0,t,x,p) (which is a C 1 
diffeomorphism) is bounded in H s ~ x . Secondly, for s sufficiently large (in the present case 
s > 5), the H s ~ x norm of the composition g o h of a mapping g of class H s ~ x and a 
diffeomorphism h of class H s ~ x can be estimated in terms of the H 3-1 norms of g and h. 
These facts follow from the results of [10] and [17]. Thus the H 3 ~ x norm of d t f can be 
estimated by an expression of the form C||/°||h s where C depends only on the if s_1 norm 
of the coefficients in the Vlasov equation. An estimate for ||T Q ^||#s-i follows. 

4. Estimates for the Einstein equations 

First the form of N ab will be examined when it is written in terms of and K ab . 
The indices of K ab will be raised and lowered using ^ ab and its inverse. 

<j) 2 L ab = V a V 6 - cfr-^VaVbcf) - t c ab V c <j) 

a2\„ .cd. 



+ 2A0«<9 t K 6c + K C b 8 t K ac ) + 4(P 2 \K ac K Cd K bd (4.1) 

t>H ab = 4>I ab + 7 cd V C 0(2V d n ab - V a K bd - V b K ad ) 

cd-r-7 ,„_ _ V7 J.„,cd\-7 „ V7 V7 J, .„C 



- V a (p-f cd V c K bd - V b (p-f cd V c K ad - V C V O 0«6 - V C V 6 <K + A 7 0K ab (4.2) 

(j)I ah = -3(R ac K c b + R bc K c a ) + 2 lab {K cd R cd ) + K abl cd R cd . (4.3) 

In (4.3) it is still necessary to make the substitution 

R ab = X(p~ 1 d t K ab + (P^VaVbcP + 2\K ac K c b + 8nT ab + AtiX^T^ - T) lab . (4.4) 

Suppose that a collection of quantities ("f a b, K ab , 0, T a/3 ) is given. These are not assumed 
to satisfy any equations. If 7a& is a solution of 

dtlab = -2<pKab, (4.5) 
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with initial datum 7 a b(0) then an estimate of the form 

\\lab{t) -7-6(0)11^. < C I \\% b (t') - lab (0)\\h + \\Vc/>(t')f Ks .4K ab (t')f Hs dt' (4-6) 

Jo 

holds for any s > 2 provided the norms appearing all exist. This can conveniently be 
proved using (A18). The time derivative of ^ a b can be estimated similarly in terms of </>, 
K a b and their time derivatives giving for any s > 3 the estimate 



Hft7-6(0llH.-i < /'||ft7a6(OllH.-i + l|V0(Olli.-i||ft««6(Oir W 

Jo 



(.1.7) 

+ \\vt(t')\\ 2 K s-i\Kb(t')\\ H s-idt'. 



Let T% c denote the Christoffel symbols of 7 a &. These satisfy 

d t T a bc = -V 6 (# c a ) - V c (#£) + V a (<j>k bc ). (4.8) 

In Sect. 5 this will be used to obtain a stronger estimate for d t t bc than could be derived 
from the estimate for d t ^ a b- Next consider the following equation which is closely related 
to (2.30). 

A f 4> = -( 7 a6 - 5 ab )d a d b - 7 ab r^V c + X[\ K \ 2 + 4n(p- 2 T 00 + 4nT]0 (4.9) 

It is shown in the appendix that if ^ a b is a metric such that \\j a b — S^Hk*- 1 + ||7a& — S a b\\ p 
is sufficiently small, if the contents of the square brackets are in L X (R 3 ) n i/ s_1 (R 3 ) and 
if A is sufficiently small then (4.9) has a solution tending to 1 at infinity. Moreover if 
bounds for \\j a b — 5 a b\\K s , ||7o6 — 8 a b\\ P and ||p||i + ||p||hs-i are given then a bound for 
Halloo + A -1 1| V</>||hs is obtained. Here p denotes the expression in square brackets in (4.9). 
Now the norms of p appearing in this estimate can be estimated in terms of the H s norm 
of K a b, the K s norm of ^ a b — S a b and the H s ~ x norm of T a/3 together with a bound on 
the size of the support of T Q/3 . If equation (4.9) is differentiated with respect to t then 
the resulting equation can be used to obtain bounds on dt4> and its spatial derivatives in 
a manner similar to the above. 

The next equation for which estimates are needed is 

-(^dtfKab + X~ 1 A 1 R ab = \-\Nab + T a b). (4.10) 

Suppose that a solution of (4.10) is given with initial data («a&(0), <9t«at>(0)). The funda- 
mental energy estimate for the equation (4.10) is obtained by multiplying it by r ) ac ^ bd d t R c d 
and integrating in space. Define 



E = 



'R3 

Then the energy estimate takes the form 



[ |0- 1 a t K a6 | 2 + A- 1 |VK o6 | 2 + |Ka6| 2 d^ 7 . (4.H) 

Jr 3 



E(t) < E(0) +c( f E(t')dt' + [ f^ 1 [r c l M (N cd + r cd )d t Rab] (t')dV, 
\Jo Jo Jn 3 



dt' 
(4.12) 
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where the constant C depends on the norm in the space L°°([0,T), AT S (R 3 )) of the quantity 
lab — 8 a b for s > 3, the norm in the space L°°([0, T), iiP(R 3 )) of the quantities n a b and 
A -1 V a 0, the L°° norm of d t 7 ab , the norm of \~ 1/2 t% c in the space L°°([0, T), fl" a_1 (R 3 )), 
the L°° norms of and and a constant A such that A _1 I < 7 < A/ where I is the identity 
matrix. This is proved in a way which is standard for quasilinear hyperbolic equations (cf. 
[16]). The important point is that the only possible A dependence of the constant C in 
(4.12) is through the quantities A _1 V</> n and A -1 ' 2 !^ and it will be seen below that the 
latter quantities can be bounded in terms of quantities which do not depend on A. 

To estimate the second term on the right hand side of (4.12) one could try to use the 
estimate 

S ab d t R ab 



J 

</R 3 



< C I \S\ 2 + \d t n\ 2 (4.13) 

R3 



for any tensor S a b- This is sufficient for N ao and for most of the terms in r a b- However it is 
not sufficient to estimate the terms involving spatial derivatives of T al3 without resulting 
in a loss of differentiability in the iteration to be carried out in Sect. 5. This difficulty can 
be overcome by the device of doing a partial integration in time already used in the proof of 
Lemma 2.1. The term to be estimated is schematically of the form /J [J R3 (<9 t KVT) (/;')] eft'. 
Now integrate by parts in time. The resulting spacetime integral can be handled by 
integration by parts in space. After partial integration in space the boundary term is 
of the form J R3 TVk. Now estimate this by the sum of a small constant times the L 2 
norm of Vk and a large constant times the L 2 norm of T a/3 . The former term can be 
absorbed into the energy and the latter can be estimated by expressing it as the integral 
of its time derivative from to t. If it is known that T% c is 0(X 1 ^ 2 ) then the estimate for 
A -1 J |V a Kb c | 2 coming from (4.12) implies a similar one for A -1 J \d a Rbc\ 2 - 

Now consider the higher derivatives of K a b- Differentiate (4.10) up to s — 1 times (using 
partial, not covariant, derivatives) and rearrange the result to give a hyperbolic equation 
for D m K, m < s — 1. Let E m be defined by replacing k in the definition of E by D m K. An 
inequality similar to (4.12) can then be obtained with E replaced by E + X)i=i an d ^ ne 
term involving and r being replaced by terms involving derivatives of N and r up to 
order s — 1. Here a bound for the H s ~ l norm of A _1 / 2 f ^ c goes into the constant appearing 
for the following reason. Many terms of the form \~ 1 D rni/ yD rn2 R with 7774,777,2 > 1 occur 
as source terms in the equation for D m K. Hence it is necessary to have bounds for spatial 
derivatives of A _1 / 2 7- These can be obtained from the equation 



Oalbc = IbdKc + IcdKb (4-14) 

Note that the quantity E + J2i=i Ei defines a norm which is stronger than ||ft a b||jy s + 
Wdt^abWH^- 1 because of the fact that it contains A -1 . 



5. Existence on a uniform time interval 

In this section the existence statement of Theorem 1.1 will be proved. First the 
definition of regular initial data must be given. 

Definition An initial data set (<7° fe (A), &° 6 (A), /°(A)) for the Vlasov-Einstein system de- 
pending on a parameter A G (0, Aq] is called regular of order s if for some p < 6 
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(i) 9ab( X ) - X6 ab belongs to the space L P (R 3 ) n K S+2 (TL S ) for each fixed A 

(ii) k® b (\) belongs to the space H S (R 3 ) for each A 

(iii) / (A) belongs to H S (TL 6 ) and has compact support for each A 

It is of course assumed that g^bW ^ s a Riemannian metric for each A. 

If a regular initial data set satisfies the conditions (i) and (ii) occurring in the hy- 
potheses of Theorem 1.1 then the quantities 7° b = X~ 1 g^ Lb and K° ab = X~ 1 k^ lb satisfy the 
conditions 

7„° 6 = Sat + OfA 1 ' 2 ), 

^ = 0(A^). (5 ' 1) 

The O-symbol is to be understood in the sense of the function spaces occurring in the above 
definition, both in (5.1) and in the hypotheses of Theorem 1.1. Consider now equation 
(2.30) on the initial hypersurface. Using the definition (3.7) the expression in the square 
bracket in (2.30) can be brought into a form which does not contain 0. The results of 
the appendix show that for A sufficiently small this equation has a solution which tends 
to 1 at infinity and that V0 is 0(A) as A — > in the space H S (H S ). It follows that if 
dtKab is defined with the help of (2.4) then it belongs to H S (R 3 ) for each A. Assumption 
(iii) of Theorem 1.1 implies that it is 0(X 1 ^ 2 ) in that space as A — > 0. In this way data 
for the reduced system (2.29)-(2.31) can be constructed from regular initial data for the 
Vlasov-Einstein system. 

To prove the existence theorem an iteration will be set up. First define 70, ko, 4>o, 
fo by extending the initial data 7 , / and the function <p° in a time independent manner 
and defining Ko(t) = k° + (dtn) t. These functions do not satisfy equations (2.29)-(2.31) 
but do satisfy the desired initial conditions. If now 7 n , K n , 4> n and f n have been defined 
the next iterate is obtained as follows. First solve (4.5) with <j> and k replaced by 4> n and 
K n and 7 as initial datum. Define 7^+1 to be equal to the solution 7. This should only be 
done on an interval [0, T n ) short enough so that 7 n +i is a Riemannian metric. Next solve 
(4.9) with 7, k, <fi and / replaced by the corresponding quantities with a subscript n and 
define n +i to be equal to the solution 0. To ensure the existence of the solution it may 
be necessary to reduce the size of T n . The value of Ao must also be restricted in a way 
described below. Now solve (4.10) with 7 replaced by 7 n and the quantities occurring in 
the definitions of N a b and r a b replaced by the corresponding quantities with subscript n. 
As initial data use k° and (c^k) I . Let K n +i be equal to the solution R obtained. Finally, 
in order to obtain / n +i solve (3.1) with 7, k and and / replaced by 7 n , K n , 4> n and f n +i 
respectively and initial datum f°. 

The next step is to show that this iteration is bounded in certain function spaces and 
that the T n can be chosen so that T n >T for all n, where T is a positive constant. Define 



a n (t) = max sup {||7 m (t) - 7m(0)|| Hs + ||3t7m(*)||jr— 1 + ||«m(*)||jy- 

+ ||ft«m(0llff.-i + \\*- 1/2 VKm(t)\\H>-i + ll/m(0llff.-i + 11^/^(^)11^-=} (5.2) 

b n (t) = max sup \\\<f>m(t)\\i> + W^Mmh + Hm(t)\\lo 
+ ||A- 1 V0 m (t)||l, s + \\X- 1/2 t m (t)\\ 2 Hs ^ + \\D s -\d t f m )(t)\\l) (5-3) 
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Lemma 5.1 Let a n (t) and b n (t) be the functions on [0,T n ) defined (5.2) and (5.3). Let 
\q be a positive constant. Define 

p n = K| 2 + 47T</>- 2 (T o)n + 4yrT n . (5.4) 

Suppose that for some p < 6 the inequalities \\^ ab — ^a&Hifs- 1 + \\-f ab — 5 ab \\ p < C\ and 
Ao(||pn(^)||oo + ||Pn(^)||H s - 1 ) < C2 froZd on the whole interval [0,T n ), where the constants 
C\ and C-i are chosen so that Lemma Al is applicable. Then there exists a constant C , 
only depending on the initial data, and positive real-valued functions D and D' on R 2 and 
R respectively which are bounded on bounded subsets such that for all n > 1: 

a n +i(t) < C + L>(supa n , sup6 n ) / a n+1 (t') + b n (t')dt' (5.5) 

Jo 

b n+1 (t) < D' (sup a n ) (5.6) 

Proof This is an application of the various estimates derived for the Vlasov and Einstein 
equations in sections 3 and 4. Note first that the required estimates for ||7 n (£) — 7n(0)||i7 s 
and ||<9t7 n (t) Wh*- 1 follow from (4.6) and (4.7) respectively. The estimates for ||ft n ||H s , 
||<9 t K n ||# s -i and \\\~ 1 / 2 \7K n \\ H s-i are obtained from (4.12) and the analogous estimate for 
E + Y,*~lEi as described towards the end of section 4. The inequalities (3.6) and (3.10) 
provide the desired estimates for ||/ n ||.ffs-i and \\dtf n \\H s - 2 - This completes the proof of 
(5.5). The estimation of <p n +i and its spatial derivatives is accomplished by applying the 
information on the solution of (4.9) given in section 4. The quantity ||A _1/,2 r n (t)||fl- s -i 
and its spatial derivatives can then be estimated using (4.8). Next apply the argument of 
the last paragraph of section 3 to estimate the last term in (5.3). Finally, the equation for 
4> m implies an estimate for it in terms of a n . □ 

I claim that the inequalities (5.5) and (5.6) imply the boundedness of a n and b n 
uniformly in n on an appropriate time interval. To see this, first let K\ be a constant 
which is greater than C and the supremum of a (t) +b (t). Next choose A and T so that 
if A and t are restricted to lie in the intervals (0, Ao] and [0, T) respectively the quantity 
||7ob — SabWK 11 - 1 + Wlab — ^ab\\ P is small enough whenever a n is less than K\ so that the 
equation for n +i can be solved and (5.5) and (5.6) hold. Let Ki be a bound for D' under 
the condition that supa n < K\. Let K3 be a bound for D under the conditions that 
supa n < Ki and sup6 n < K 2 . Reduce the size of T if necessary so that 

(C + K 2 K 3 T)e KaT < K x . (5.7) 

By induction a n (t) < K\ and b n (t) < K 2 for all t G [0, T) and all n. Thus the iteration is 
bounded as claimed. This bounded iteration is what is needed to show the existence of a 
regular solution of the equations, as will now be shown. 

Definition A solution (<7afc(A), fc a &(A), 0(A), /(A)) of the Vlasov-Einstein system is called 
regular of order s if for each A 

(i) g ab {\) - \5 ab and its time derivative belong to L°°([0, T), if S (R 3 )). 

(ii) k ab (X) belongs to L°°([0, T), H S (K 3 )) and d t k ab (X) belongs to L°°([0, T), i7 s_1 (R 3 )) 
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(iii) 0(A) and its time derivative belong to L°°([0, T), if S+1 (R 3 )) 

(iv) /(A) and its time derivative belong to L°°([0, T), i7 s_1 (R 6 )) 

The boundedness of the iteration shows that 7™ is bounded in C 1 ([0, T), if s (R 3 )) 
for each fixed A. In particular 7^ is bounded in L°°([0, T) x R 3 ) and V7 n is bounded 
in L°°([0, T), ii" s_1 (R 3 )). Each of these spaces is the dual of a Banach space and so any- 
bounded sequence has a subsequence which converges in the weak* topology. By passing to 
a subsequence again it can be ensured that the time derivatives of the 7™ converge in that 
sense. In the following we will always use the same notation for a subsequence as for the 
original sequence. Now let R > be a real number and let 7 n (-R) denote the restriction of 
7 n to the closed ball Br of radius R about the origin in R 3 . A similar notation will be used 
for other objects defined on R 3 . The sequence {■~f n (R)} is bounded in C 1 ([0, T), H s (Br)). 
Furthermore the embedding of H s (Br) into H s ~ 1 (Br) is compact. Using the vector- 
valued Ascoli theorem[9] it can be concluded that there is a subsequence which converges 
strongly in C°([0, T), H s ~ 1 {Br)). A diagonal argument shows that the a subsequence can 
even be chosen so that this is true for all R. By passing repeatedly to a subsequence and 
using the same argument it is possible to obtain similar convergence statements for k n , 
V(/> n - For f n it is even possible to get strong convergence in C°([0, T), if s_2 (R 6 )) since the 
supports of all f n are known to be contained in a fixed compact subset of R 6 . In the end 
we obtain a sequence which converges weakly in a space defined globally on R 3 (or R 6 in 
the case of f n ) and strongly in spaces defined on compact subsets. The weak convergence 
ensures that the limit of the this sequence belongs to the spaces specified in the definition 
of a regular solution while the strong convergence ensures that it does define a solution of 
the Vlasov-Einstein system. 

6. Convergence to the Newtonian limit 

First it is necessary to get hold of the functions Jn and U which occur in the statement 
of Theorem 1.1. Now /(A) is bounded in C 1 ([0, T), iP _1 (R 6 )). Choose a sequence {A n } 
converging to zero. Let f n = f(X n ). (This should not be confused with the sequence of 
iterates used in the last section.) The sequence {f n } is bounded in C 1 ([0, T), i7 s_1 (R 6 )) 
and hence, the Ascoli theorem for vector-valued functions, has a subsequence which con- 
verges strongly in C°([0, T), i7 s_2 (R 6 )). Call the limit Jn- Next consider the sequence 
{A~ 1 (l — 4> n )} where 4> n = <fi(X n ). This is bounded in C 1 ([0,T) x R 3 )) and the sequence 
{A~ 1 V0 n } is bounded in C 1 ([0, T), H S (R 3 )). By the type of argument used in the previ- 
ous section it can be seen that, after passing to a subsequence, {A~ 1 (l — 4> n )} converges 
uniformly on compact subsets to a limit, which will be denoted by U, and that for any 
R > the restriction of A~ 1 V0 n to Br converges in C°([0, T), H s ~ 1 (Br)) to the restric- 
tion of VU . Moreover, using the weak* convergence argument, U is bounded and VU is 
in L°°([0, T), H S (R 3 )). In the same way it can be concluded that there is a sequence {A n } 
converging to zero such that 7 n = 7(A n ) and K n = K(X n ) converge in suitable function 
spaces as n — > 00. From equation (4.4) it can be seen that the Ricci tensor of 7™ tends 
to zero as n — > 00. Hence the limiting metric is flat. Consider now the time derivative of 
the connection. This is given by equation (4.8). We know already that V0 is O(X) and 
that Vk is 0(X 1 / 2 ) in certain spaces. Hence equation (4.8) implies that T£ c is 0(X 1 / 2 ). 
Integrating in time shows that the Christoffel symbols of ^ a b are 0{X 1 / 2 ). Now the partial 
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derivatives of ^ ab with respect to the spatial coordinates can be written in terms of these 
Christoffel symbols and 7 a 5 itself as in (4.14) and it follows that the partial derivatives 
of must be 0(X 1 ^ 2 ). In particular the partial derivatives of the limit of the 7 n vanish 
and so this limiting metric must in fact be given by S ab . As a consequence the limit of 
the K n must also vanish. The characteristics of the Vlasov equation converge uniformly 
to those of the non-relativistic Vlasov equation with force term VU along the sequence 
A n . Hence /n must coincide with the unique solution of the latter equation with the C 1 
initial datum /j^. Passing to the limit in (2.30) and (3.1) then shows that and U 
satisfy the Vlasov-Poisson system with initial datum Unfortunately the convergence 
statements derived up to now are not enough to ensure that the function U satisfies the 
standard boundary condition that U — ■> as r — > oo. On the other hand it is known that 
U is bounded so that it can be made to satisfy the boundary condition by adding to it 
an appropriate constant. It will be supposed from now on that this alteration has been 
made so that U vanishes at infinity. Since the solution of the Vlasov-Poisson system with 
a given C 1 initial datum is unique, it follows that for every sequence {A n } the sequence 
{A _1 V</> n , f n } has a subsequence tending to the same limit, namely (VU, fx). Hence the 
restriction of A _1 V<^ to any ball of radius R converges to V£7 in H s ~ l and f n converges 
to / in H s ~ 2 . It is also possible to show weak* convergence globally in R 3 . 

The convergence statements obtained up to now can be improved by estimating the 
difference between the Newtonian and relativistic solutions. 

a / (a- 1 (^ - 1) + u) = -( 7 ab - 5 ah )d a d h {\-^) - 7 a& r: b (A- 1 a c ^) 

+ [\k\ 2 + 47r</r 2 Too " 47rT 00 (0) + 4ttT]0 



Hence 



|| A 1 {4> - 1) + 17 < C(||7a6 - S ab \\ K s-i + ||K a &||#.-2 
+ ||Too-Too(0)|| Hs -2+A) 

Differentiating (6.1) gives the estimate 



(6.2) 



rs-l 



||A <p + U\\ K ' < C(\\"/ ab - S ab \\ K s-i + \\dt<yab\\H° 

+ \\ K ab\\H s ~ 2 + \\9tK a b \\h 3 ~ 2 + 11^00 ~ Tq \\ H s-2 + A) 



(6.3) 



Using the definition of T 00 it can be shown straightforwardly that the terms involving T Q0 
in (6.2) and (6.3) can be replaced in these inequalities by expressions involving f — /n- 

||A -1 ((/> - 1) + U\\ K ° < C(||7a6 - 5 ab \\ Ks -i + \\K ab \\ H s-2 
+ Wf-fN\\ H ->+X) 

The inequality (6.3) can be modified similarly. Subtracting the Vlasov equation for 
from that for / gives the estimate 

ll/W " fN(t)\\ 2 H s- 2 < ||/(0) - M0)\\h.-» +CJ\W) - fN(t')\\ 2 Hs - 2 

+ habit') -S ab \\ 2 K3 -2 + \\\- 1 v c <f>(t') + V c ?7(OIIh-i +X 2 )dt' 
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and an analogous estimate for \\f(t) — /at(*) || jy«- 2 • Equation (2.29) implies an estimate of 
the form 



hab(t) -7a6(0)||| f8 -i < C I \\ lab {t') -7a6(0)||? f8 -i + K&(0 llff.-i^' (6-6) 

Jo 

To close the argument and obtain a useful differential inequality, it remains to estimate K ab . 
This can be done by examining carefully the third term in (4.12). This is mostly routine 
but there are two expressions which require particular care and these will be handled 
explicitly here. The first, which arises from N ab is: 



l„,ac„,6d. 



A 7 7 



l V C V d <j)d t K ah dV 1 



R 3 



< 



l„,ac„,bd , 



A 7 7 



R3 



i vJv c {d t K ab )dv 1 



+ ... 



< CA^diA-iv^Hl + HA-^v.O^Ui) + . . . 

The other, which arises from r ab , is | J R3 7 ac 7 6d </> _3 Too7cd<9^ a ; ) <iV r 7 | To estimate this, first 
note that T 00 — T 00 (0) can be controlled straightforwardly so that T 00 can be replaced 
without loss of generality by T oo (0). Now the fact can be used that t 00 (0) = -d a T Oa (0). 
After this substitution has been made it suffices to integrate by parts in space. 
The quantity dt^ ab satisfies U^V/abHi?*- 1 < H^atH/f^- 1 - Furthermore: 



hab(t) ~ S ab \\ K s-i < ||7o6(*) -7o6(0)||h—i + ||7o6(0) -^Hh— i 



(6.7) 



Define 

a(t) = SUp{||7 ab (t) - 7 a6 (0)||^ s _i + ||Ka6(*)||H— i + ||&«a&(*)||ff.-2 
A 

+ WX-^VKatmH-' + ll/(*) " fN(t)\\%-' + W d tf(t) ~ dtfN(t)\\ 2 H s-i} 

Then the above estimates show that 



(6.8) 



a(t) < C(A 1/2 + / a{t')dt') 



(6.9) 



It follows using Gronwall's inequality that a(t) = 0(X 1 ^ 2 ). 

The meaning of the order symbols in the statement of Theorem 1 can now be explained. 
They refer to the function spaces obtained from those occurring in the definition of a regular 
solution of the Vlasov-Einstein system by replacing s by s — 1. In other words, convergence 
is obtained in a space involving one less derivative than that where the existence of the 
solution has been obtained. Given this definition, the conclusions (i)-(vi) of Theorem 1.1 
follow from (6.9), (6.7) and (6.4). 



7. Solution of the constraints 
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In terms of 7 a & and n ab these take the form 

R 7 _ \\ K \ 2 = 16ttA(/)- 2 Too, (7.1) 
7 a& V a « 6c = -8^-^00 (7.2) 

for a maximal hypersurface. Let fi = 4>~ 2 T 00 and J a = (j)~ 1 T 0a . Note that when fi and J a 
are expressed in terms of / and 7 a & using (3.7) then it is seen that they do not explicitly 
depend on <fi. To solve these equations we start with the following A-dependent objects: 

(i) a function jj of the form /2(0) + 0(X 1 ^ 2 ) 

(ii) a vector J a which is 0(A) as A — > 

(iii) a metric ^ ab satisfying ^ ab = 5 ab + 0(A 3/2 ) 

(iv) a traceless symmetric tensor R ab satisfying R ab = 0(X 1 ^ 2 ) and 

It is assumed that fi G _£P(R 3 ), J a G if s (R 3 ) and 7 a 6 — G if|(R 3 ) for some 5 in the 
interval (—1, —1/2) and that all these objects depend continuously on A in the given spaces. 
The O-symbols here also refer to those spaces. Furthermore, it is assumed that ft and J a 
are compactly supported. Tensors k ab of the type required in (iv) can be contructed in a 
standard way using the York decomposition. From these objects it is possible to determine 
initial data for the Einstein equations by solving the equation [3] 

A^V + (l/8)(-i^V + A|k|V" 7 ) + 27rA/}-i/>" 3 = (7.3) 

and definining 7 = -0 4 ^, k = i/j~ 2 k, 4>~ 1 T 0a = J a = ifj~ 6 J a and 4>~ 2 T 00 = fi = i/j~ 8 fl. It is 
well known that equation (7.3) can be solved for a unique ip which tends to 1 at infinity 
and which has the property that ifj — 1 belongs to a weighted Sobolev space provided that 
Tab is close to 5 ab in a weighted Sobolev space. Because of assumption (iii) above this 
will be satisfied when A is close to zero. Note that the condition 5 > — 1 implies that if 
s > 2 the space i/|(R 3 ) is continuously embedded in L P (R 3 ) for some p < 6. Let ipi 
be the solution of the equation Afipi = —2rcjl(0) which tends to zero at infinity. Then 
since U satisfies the Poisson equation and ip is identically one when A = it follows that 
ipi = (1/2)U. Now a comparison with (7.3) shows that 

A / (V-A^i)=0(A 3 / 2 ) (7.4) 

Hence V = 1 + A^i+0(A 3 / 2 ). This can be used to calculate the Ricci tensor of 7^ to order 
A using the formula for conformal transformations. The result is that the contributions of 
order A on the right hand side of (2.4) cancel so that assumption (iii) of Theorem 1.1 is 
satisfied by the given initial data. Now the significance of this assumption can be discussed. 
In [21] it was shown that a necessary condition for the existence of a regular Newtonian limit 
is that the coefficient of A 2 in the expansion of the spatial part of the metric, considered 
as a linearised metric, has vanishing linearised Bach tensor. This means that it satisfies 
the linearised version of the condition of conformal flatness. Another way of expressing 
this would be to say that the spatial metric is conformally flat to second order in A. It 
was in order to ensure that this condition was satisfied that 7 was chosen to be flat up to 
first order in A. The computation which has just been done shows that the role of this 
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condition in the present paper is to ensure the correct behaviour of d t n a b for A — > on the 
initial hypersurface in order to make the iteration work. 

To get a complete initial data set for the reduced equations it is still necessary to 
describe how the initial values for / are to be obtained. This is done as follows. Suppose 
that a test metric 7 & has been chosen as in (iii) above. Let / be a non-negative real- valued 
function of compact support on the mass shell defined by 7 a &. Define \x and J a in terms of 
/ and 'jab in the same way as \i = <p~ 2 Too and J a = </> _1 To a are defined in terms of / and 
•jab (i.e by (3.7)). Now construct the quantities 7 a 6, K a b, To a and Too as described earlier 
in this section. Let p a = if; 2 p a and define f(p a ) = ifj~ 8 f(p a ). It then follows that Too and 
T 0a arise from / via the equation (3.7) and a full initial data set for the Vlasov-Einstein 
system is obtained. 



Appendix 1. Some elliptic theory 

The discussion here will be limited to results which are not standard. For background 
material on the Poisson integral see [7] . In this appendix A always denotes the Laplacian 
of the standard flat metric on R 3 . For any function on R 3 and a positive real number e 
let <j> e (x) = 4>(ex). If / is a continuous function belonging to L P (R 3 ) n L°°(R 3 ) for some 
p < 3/2 then it possesses a Newtonian potential u which is C . This is a solution of the 
Poisson equation Au = f in the sense of distributions and can be obtained as the limit for 
e — > of the Newtonian potentials of the compactly supported functions 4> e f where <fi is 
any C°° function of compact support which takes the value 1 in a neighbourhood of the 
origin. The Newtonian potentials of the functions <p e f can be represented by the familiar 
Poisson integral. They tend to zero at infinity. Pick a sequence e n tending to zero and let 
fn = 4>£ n f- Then the f n are bounded in L°°(R 3 ) and converge to / in L P (R 3 ). Let u n 
denote the Newtonian potential of f n . A useful estimate for the Poisson integral will now 
be derived (cf. [2]). For any R > 0: 

R3 F — y\ J\x-y\<R \ X ~ U\ J\x-y\>R \ X ~ V\ 

^cdi/iu^ + n/ii^ 3 -^). 

Here it has been assumed that the conjugate exponent q of p is greater than 3, which 
implies that p < 3/2. Putting R = (||/|| P /||/||oo) p ^ 3 gives the estimate 

lkn||oo<C||/ n ||^/ 3 ||/ n ||^/ 3 , p<3/2. (Al) 

The derivatives of u n can be estimated similarly. 

\\Vu n U<C\\f n \\l^f n \\^'\ p<3. (A2) 

These estimates show that u n and its first derivatives converge uniformly to u and its first 
derivatives. In particular u tends to zero as \x\ — > 00. Now all the functions u n have the 
property that u n (x) = 0(|x| _1 ) and V«„(x) = 0(|x| -2 ) as \x\ — > 00. Hence a partial 
integration shows that 

ll^^nl^ — ||^n||oo WfnWl- (A3) 
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If / e L 1 (R 3 ) it follows that V«„ is a bounded sequence in L 2 (R S ). Since this sequence 
also converges pointwise it follows from Fatou's lemma that V« is in L 2 (R 3 ) and that it 
satisfies the equivalent of (A3). Putting this together with (Al) gives 



V«|| 2 < CH/H^ 



11/11? 



5/6 



(A4) 



Another estimate which is satisfied by u when / e H S (H S ) is 



\\d a d b u\\ H s < \\f\\ H s. 



(A5) 



It has now been shown that if / e L 1 (R 3 ) flif s (R 3 ) then its Newtonian potential u belongs 
to L°°(R 3 ) and d a u belongs to if s+1 (R 3 ). In the following another related result will be 
required. Let H ab be a tensor on R 3 whose components belong to L P (R 3 ) n L°°(R 3 ) 
for some p < 6 and which has the property that d c H ab belongs to i7 s_1 (R 3 ). Let h 
be a function which is bounded and whose first derivatives are in # S (R 3 ). Consider the 
equation Au = H ab d a dbh. By Holder's inequality the expression on the right hand side 
of this equation belongs to L g (R 3 ) for some q < 3/2. Thus it can be concluded from the 
above discussion that u is in L°°(R 3 ) and an estimate for its norm in that space follows 
from (Al). Now let Bf = (j) en H ab and let u n be the solution of Au n = Hfd a d b h. 



Thus the L 2 norm of Vtt n is bounded and it is possible to argue as above that Vw is in 
L 2 (R 3 ) and that its norm can be estimated in terms of || 1| oo , ||<9a-ff ab ||2 and ||<9b/i||2- 
Lemma Al Let 'jab be a Riemannian metric such that 'jab — $ab G L P (R 3 ) n if s (R 3 ) for 
some p < 6 and some s > 3 and let F a be its contracted Christoffel symbols. Let p be a 
function belonging to L 1 (R 3 ) n if s_1 (R 3 ) and A a non-negative real number. Then there 
exist positive constants C\ and Ci such that if \^ a b ~ ^ab\\ P + ||7o6 — ^abWK 8 - 1 < C\ an d 
(IIpIIi + IIpIIh 8 - 1 )-^ — ^2 then the equation 



has a unique solution with the property that <fi — > 1 as \x\ — > oo. Moreover for any k < s 
the solution satisfies an estimate of the form 




= (5 1 



ab 



(A6) 



A _1 (0 - l)||oo + ||A~ 1 V0||i 3 -fc <C 



(A7) 



for a constant C only depending on C\, and \\iab ~ £>ab\\K a - 
Proof Define an iteration by solving 



A0 n+ i = (5' 



ab 



7 a6 )a a a 6 0n + r a a a n + Ap0, 



(A8) 



•n 
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with <j) n+ i — > 1 as \x\ — > oo and O = 1. Let g n = T a d a (p + Xp4>. If 4> n £ L°°(R 3 ) and 
V0 n G _£P(R 3 ) then g n e L 1 (R 3 ) fl if s_1 (R 3 ) and so by the above remarks concerning 
the Poisson equation the solution n +i exists. Furthermore 

||?n||l + Iknllffk-i < C(ll7a6 ~ S ab \\ K k + A(||p||i + ||p|| H *-i )) ( ||0 n || ^ + ||V0 n || H fc) (A9) 

for any k with 3 < k < s. Thus there exist constants C\ and C 2 such that if the inequalities 
||7o6 - 5 ab \\ LP + W'jab ~ SabWics < C\ and (||p||i + ||p|| H s-i)A < C 2 hold then 

||0n+l 4 > n || 00 + ||V0 n+ i - V0 n ||H3 < K(\\(j) n - </> n -i||oo + ||V^ n - V(j) n -i\\ H 3) (A10) 

for some K < 1. It follows that {0 n } and {V0 n } are Cauchy sequences in L°°(R 3 ) and 
if 3 (R 3 ) respectively and hence n converges to a solution of (A6) with — > 1 as |x| — > 00. 
This solution satisfies 

Woo + ||V0|| H 3 <C (All) 
for some C depending only on C\ and C 2 . Dividing (A6) by A gives 

A(A- V) = ($ ab - r^dad^x- 1 ^) + r a d a (\- V) + P</> (A12) 

Estimating the quantity \~ 1 (p n in the same way as <J) n was estimated above shows that 
{X~ 1 ((/) n — 1)} and {A -1 V</> n } are Cauchy sequences in L°°(R 3 ) and if 3 (R 3 ) respectively 
and that (A7) holds for k = 3. To see the that it is true for any k < s, use the information 
we already have in the right hand side of (A12). Using the estimates for the Poisson 
equation stated at the beginning of this appendix then shows that if the hypotheses of the 
lemma hold and if s > 4 we obtain (A7) for k = 4 after possibly reducing the size of C\. 
This process can be repeated until k = s. □ 



Appendix 2. Estimates for modified Sobolev spaces 

This appendix is concerned with proving some useful estimates for the modified 
Sobolev spaces K s (R 3 ) introduced in section 3. Recall that 

ll/lk. = ll/lloo + ||V/|| H .-i. (A13) 

The results to be proved are analogues of the results (3.2) and (3.3) for functions belonging 
to ordinary Sobolev spaces and are proved in a similar way. Estimates of this type are 
discussed in [14], [16] and [23]. First analogues of (3.2) will be discussed. Suppose that /, g 
belong to K S (R 3 ) for some s > 2. Note first the obvious fact that ||/y||oo < ||/||oo||<7||oo- 
If a is a multi-index with 1 < \a\ < s then 

D«(fg) = D a fg + fD a g + £ -^D^fD^g. (A14) 

The sum is taken over all multi-indices (3 and 7 with \(3\ + (7] = \a\ and max(|/3|, I7I) < \a\. 
The first estimate we wish to prove is that 

\\D a (fg)h<C\\f\\ K .\\9\\K. (A15) 
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for some constant C. It is elementary to estimate the first two terms in (A14) by the 
right hand side of (A15) and so we can concentrate on the third term. Consider one of 
the summands there. Suppose first that either (3 or 7 is less than s — 3/2. Without loss of 
generality we can assume that it is (3. Then 

WD^fD-ygh < H^/IUlZ^lh < CII^/II^HD^IIa (A16) 

If on the other hand both (3 and 7 are greater than s — 3/2 then it can be concluded that 
s < 3. The only estimate which remains to be done to establish (A15) is 

\\DfDg\\ 2 < \\DfU\Dg\U < \\Df\\ m \\Dg\\ m , (A17) 

where the Holder and Sobolev inequalities have been used. The inequality (A15) shows 
that multiplication is a continuous mapping from if s (R 3 ) to itself. This is a statement 
of a weaker type than (3.2). A stronger result could presumably be obtained using the 
Gagliardo-Nirenberg inequality but that will not be attempted here since (A15) is sufficient 
for the applications in this paper. In a similar way it can be shown that if / e fP(R 3 ) 
and g e X S (R 3 ) for s > 2 then fg G H S (K 3 ) and 

\\fg\\H. < \\f\\H>\\g\\K.. (Ais) 

An estimate related to (3.3) can also be obtained if s > 3, namely 

\\D a (fg) - fD a g\\ 2 < C\\f\\ K >\\9\\K'-i, 1 < \a\ < s (A19) 

The final estimate which is required concerns composition with a C r function. Suppose 
then that U is an open interval in R and F : U — > R a C r function whose derivatives up 
to order r are bounded on U. Let / be a function in K S (R 3 ), where 2 < s < r, whose 
range is contained in U. Now 

D a (F(f)) = J2Cra 1 ...a l ^D^f...D^f (A20) 

where 1 < k < \a\, \a± \ + . . . + \a\\ = \a\ and \oci\ > 1 for all i. By an argument similar to 
those given above it is seen that the terms on the right hand side of (A20) can be estimated 
straightforwardly unless s = 2. In that exceptional case the embedding i7 1 (R 3 ) — > L 4 (R 3 ) 
can be used. There results the estimate 

\\F(f)\\ K s<C\\F\\ C r\\f\\ s K3 . (A21) 

In this paper (A21) is only needed to estimate (/> _1 and the inverse metric r ) ab . In both 
cases it is applied to the function F(f) = 1/ / with U = (c, 00) for some c > 0. 
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